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A detailed theoretical analysis of the Gaussian fluctuations of the order parameter in 
layered superconductors is performed within the Ginzburg-Landau (GL) theory. The 
available results for the Gaussian fluctuations are systematized and a number of novel 
formulae for the fluctuation magnetization, nonlinear magnetic susceptibility, heat ca- 
pacity and high-frequency conductivity in layered superconductors are derived. We 
propose several new prescriptions: how to determine the life-time constant of fluctu- 
ation Cooper pairs tq, the in-plane coherence length i;ab{0), the energy cutoff parameter 
= ch? /2mabiab(^)' ^'^'^ Ginzburg number Eqj. It is demonstrated, for example, 
how the spectroscopy of the life-time can be used to verify the existence of depairing 
mechanisms in layered cuprates. The ultraviolet regularization of the GL free energy is 
carried out by means of the well-known from the field theory (^-function method. We 
further show that the archetype of the latter method has its origin in the century of en- 
lightenment and the novel result is that the fiuctuation part of the thermodynamic vari- 
ables of the layered superconductors can be expressed in terms of the Euler F-function 
and its derivatives. Universal scaling curves for the magnetic field dependence of the 
paraconductivity aab{Tc, B), fluctuation magnetization M{Tc,B) and the heat capacity 
C{Tc,B) are found for the quasi-2D superconductors at Tc and further related to the 
form-factor of the Cooper pairs. 

Keywords: Gaussian fluctuations, Ginzburg-Landau theory, layered superconductors 

1. Introduction 

The writing of this review is provoked by the progress in the study of fluctuation 
phenomena in high-Tc superconductors.^ The small coherence lengths in the layered 
cuprates £,a — S,b give rise to a very high density of fluctuation degrees of 

freedom oc l/(^ab('-')Cc(0)) which makes the fluctuation effects easier to be observe 
in the high-Tc- rather than in the conventional superconductors. An intriguing 
feature of the fluctuation effects to be pointed out is that they can be observed 
even in the case when the interaction between fluctuations is vanishing or can be 
treated in a self-consistent manner. In such a case, for high quality crystals, the 
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fluctuations are of Gaussian nature and their theory is very simple. A number of 
good experimental studies have already been performed in the Gaussian regime thus 
initiating the Gaussian fluctuation spectroscopy for high-Tc materials. 

By spectroscopy here we imply only those experiments with trivial theory where 
every measurement provides an immediate information for some parameter(s) im- 
portant for the material science or fundamental physics of these interesting mate- 
rials. Half a century ago Landau used to speak about himself as being the greatest 
trivializator in the theoretical physics. At present, the Ginzburg-Landau (GL) 
theory (called by Ginzburg also ^f-theory) is the adequate tool to describe the fluc- 
tuation phenomena in the superconductors. The parameters of the ^'-theory, such 
as coherence lengths, relaxation time tq.iij of the GL order parameter the GL 
parameter k^^ = Xab{0) / £.ab{Q) are also "meeting point" between the theory and 
the experiment. 

Prom one hand, these parameters are necessary for the description of the ex- 
perimental data and from another hand they can be derived form the microscopic 
theory using the methods of the statistical mechanics. That is why the determi- 
nation of the GL parameters is an important part of the investigations of every 
superconductor and the Gaussian fluctuation spectroscopy is an indispensable tool 
in these comprehensive investigations. 

The purpose of this review is to systematize the known classical results for the 
GL Gaussian fluctuations, to derive new ones when needed, and to finally give 
suitable for coding formulae necc^ssary for the further development of the Gaussian 
spectroscopy. The derivation of all results is described in detail and trivialized to the 
level of the Landau-Lifshitz encyclopedia on theoretical physics,^ the textbooks by 
Abrikosov^ and Tinkham^ or the well-known reviews by Cyrot'' on the GL theory, 
by Bulacvskii'' concerning the layered superconductors with Josephson coupling, 
and by Skocpol and Tinkham^ on the fluctuation phenomena in superconductors. 
The present work is intended as a review on the theoretical results which can be 
used by the Gaussian spectroscopy of fluctuations but no historical survey of the 
experimental research is attempted. Therefore we do explicitly refer to only a 
limited number of experimental studies in this fleld. Instead, the reader is referred 
to the citations-reach conference proceedings,^ but even therein a number of good 
works are probably not included. We do not cite directly even the epoch-creating 
paper by Bednorz and Miiller but its spirit can be traced to every contemporary 
paper on high-Tc superconductivity. Even to focus on the theoretical results related 
to fluctuation phenomena is a very difficult problem by itself and therefore, when 
referring to any result one should imply "to the best of our knowledge...". One of 
our goals was also to fill the gap between the textbooks and experimentalists' needs 
for a compilation of theoretical formulae written in common notations, appropriate 
for direct use. 

Of course, there is a great number of interesting physical situations especially 
related to vortices where the fiuctuations are definitely non-Gaussian. Those prob- 
lems fall beyond the scope of the present review and we include only some references 
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from this broad field in the physics of superconductivity.^ 

The review is organized as follows: in Sec. 2 the case of weak magnetic fields 
is considered and the thermodynamic variables are expanded in power series in 
the dimensionless magnetic field h = B^/ Bc2{0)- The standard notations for the 
thermodynamic variables in a layered superconductor arc then introduced in Subscc. 
2.1, and Subsec. 2.2 is dedicated to the Euler-MacLaurin summation formula in the 
form appropriate for the analysis of the GL results for the free energy and its 
ultraviolet (UV) rcgularization. A systematic procedure to derive the results for a 
layered superconductor from the results for a two-dimensional (2D) superconductor 
is developed in Subsec. 2.3 and the action of the introduced "layering" operator 
L is illustrated on the example of the formulae for the paraconductivity. Further 
we consider the static paraconductivity in case of perpendicular magnetic field as 
well as the high frequency conductivity in zero magnetic field. The power series for 
the nonlinear magnetic susceptibility and the magnetic moment in the Lawrence- 
Doniach (LD) model arc derived in Subscc. 2.4 and the e-mcthod for summation 
of such divergent series is described in Subsec. 2.5. For practical purposes a simple 
FORTRAN90 program is given in the Appendix. 

Further in Subsec. 2.6 we present the power series for the differential suscep- 
tibility and general weak-magnetic-field expansion formulae for the magnetization. 
Section 3 is dedicated to the study of the strong magnetic fields limit. Firstly, in 
Subsec. 3.1 the general formula for the Gibbs free energy in perpendicular to the 
layers magnetic field is analyzed. The fluctuation part of the thermodynamic vari- 
ables is found then by differentiation in Subsec 3.2. Subsec. 3.3 is devoted to the 
self-consistent mean-field treatment of the fluctuation interactions in the LD model. 
The important limit case of an anisotropic 3D GL model is considered in Subsec. 3.4 
where we derive the Gibbs free energy and the fluctuation magnetic moment. In 
Sec. 4 an account is given of the fitting procedure for the GL parameters which 
rests on theoretical results and some recommendations for the most appropriate 
formulae are also given for determination of the cutoff energy in Subsec. 4.1, the 
in-plane coherence length ^a6(0) in Subsec. 4.2, the Cooper pair life-time constant 
To in Subsec 4.3, and the 2D Ginzburg number in Subsec. 4.4. All new results 
derived throughout this review are summarized in Sec. 5 and some perspectives for 
the Gaussian spectroscopy are discussed as well. 

2. Weak magnetic fields 
2.1. Formalism 

Before embarking on a detailed analysis we shall briefiy introduce all entities en- 
tering the basic for our further considerations quantity — the GL functional G for 
the Gibbs free energy in external magnetic field H^'"'*-'. For compliance with the 
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Fig. 1. Pictorial representation of material parameters that enter Eq. (1). The z-coordinates of 
the Cu02 planes, satisfying periodic boundary conditions, are 2jv+i,n = ^i,n+i, where n labels 
the unit cell and j = 1, . . . , A'^ is the index of the Cu02 plane within each unit cell. Periodicity 
in c-direction is designated by s, i.e. < Zj^ < s. ^ab{0) stands for the in-plane coherence 
length extrapolated to T = 0. The Josephson coupling energies between the neighboring planes 
Jj = aQ-fj are parameterized via dimensionless quantities -jj. Lastly, the effective charge and the 
in-plane effective mass of the Cooper pairs are, respectively, |e*| = 2|e| and m„(,. 



previous works we follow the standard notations in which G reads 
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with A being the vector potential of the magnetic field B = V x A. 

The material parameters in this sizable expression are illustrated in Fig. 1, thus 
we only need to note that the GL potential a{e) = age is parameterized by Oq = 
h'^/2mabCbi^): and e = ln(r/rc) w (T - Tc)/Tc is the reduced (dimensionless) 
temperature. If not otherwise stated we shall make use of the SI units, thus the 
magnetic permeability of vacuum fxo — Att x 10"''. 

Here we will restrict ourselves to the study of fluctuations in the Gaussian regime 
in the normal phase not too close to the critical line Hc2 (T) . In this case the nonlin- 
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ear term in G[^, A] is negligible, 6|\l/|'*/2 — + 0. For the normal phase the magnetiza- 
tion is also very small and with high accuracy iiqH.^'^^^^ B = ;Uo(H + M) « /ioH. 

To begin with, consider the simplest case of zero external magnetic field H^"'^*^ = 
0. Given the above assumption for "if, the GL functional is a quadratic form and 
one needs to sum over all eigenvalues of the energy spectrum 



2mab 



+ £cj{Pz), £cj{Pz) 



(2) 



where ScjiPz) are the tight-binding energy bands describing the motion of Cooper 
pair in z (c) direction, p = {px,py) is the in-plane (afe-plane) momentum of the 
fluctuating Cooper pairs and 9 = Pzs/h G (0, 27r) is the Josephson phase. For 
a single layered material, N = 1, this corresponds to the well known Lawrence- 
Doniach model, ^ 



a;(L°)(6i) = 27i(l-cos6i) 



(3) 



while the case A'' = 2 is the Maki-Thompson (MT) model, proposed independently 
by Hikami and Larkin^^ as well. 



(MT) 



(e) = 71 + 72 + {-ly^Jll + 72 + 27172 cos^, j = l,2. 



Thus, the sum over the energy spectrum gives 



(4) 



(5) 



where ^p,p^,j is the wave function of the superconducting condensate in momentum 
space. We use standard periodic boundary conditions for a bulk domain of volume 
V = Lx X Ly X Lz which give 
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In order to calculate the fluctuation part of the Gibbs free energy G'(r) at zero 
magnetic field, one usually solves for every point in the momentum space p,Pz,j 
the Gaussian integral 
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where * = = + i'i", and if G (0, 2it). Making use of this auxihary result 

the calculation of the fluctuation part of the Gibbs free energy reduces to summation 
over the spectrum of an effective Hamiltonian, i.e. 



G 
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(8) 
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or, taking into account Eqs. (2) and (6), 

{p'^/2mab) + aouj'-'^\0) + ao€ ao 1 

• (9) 



G , ^ f d (7rp2) 1 ^ / rf6) 

.7=1 



In view of the further calculations it is also useful to introduce a dimensionless 
in-plane kinetic energy 

x=i: = ^ e(0,c), (10) 



2mabao \ h 

bound by a dimensionless cutoff parameter c which we consider to be an important 

parameter of the GL theory when appHcd to copper oxide superconductors. Later 
in Sec 4.1 we demonstrate how the value of the dimensional cutoff energy e^^, 

can be determined by fitting to the experimental data. An immediate simphfica- 
tion to Eq. (9) can be achieved by dropping the ^k^T/ao const multiplier in 
the argument of the logarithm as it is irrelevant for the critical behavior of the 
material. Furthermore, since fluctuational observables are related to non analytical 
dependence of the Gibbs free energy on the reduced temperature, we can substitute 
T = Tc(l + e) Tc and the free energy per unit volume F{e) is cast in more elegant 
form, 

F(e)^ = — Ir-rr— / dx—y(b—\n(x + uj) '(9)+e), (12) 

L^LyL, 4<2^(0) s io Njr[J^T^ ^ y> )^ y > 

that could easily include the (1 + e)-factor in all cases when necessity appears. The 
physical meaning of this important for our further considerations expression is fairly 
transparent: one has to integrate with respect to the Josephson phase 9, which 
describes the motion of Cooper pairs in c-direction, and to take into account as 
many different Cooper pair energy bands as are there the different siiperconducting 
layers per unit cell. Finally, integration with respect to the in-plane Cooper pair 
kinetic energy is to be carried out. 

Consider now the important case of an external magnetic field applied parallel 
to the c-direction, i.e. perpendicular to the Cu02 planes, B = (0,0,5). In this 
case, the in- plane kinetic energy of the Cooper pairs acquires oscillator spectrum, 
corresponding to the quantum mechanical problem of an electron in an external 
magnetic field, ^ 

^ Swc , (13) 



2mab V 2 
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where n = 0,1, 2, 3,... is a non- negative integer and ojc = |e*| B/rriab is the cy- 
clotron frequency. The integration over the momentum space is thus reduced to 
summation over oscillator energy levels 



" (14) 



n=0 



where ric = c/2h and $o = 27r?i/ |e*| = 2.07 fTm^ is the flux quantum. The energy 
cutoff is to be applied now to the oscillator levels,^ hudric + 5) = ccq. Let us recall 
that the equation for the upper critical field Hc2 (T) within the GL theory is nothing 
but the equation for annulment of the lowest energy level, ^/iWc + a(e) = 0. Thereby 
introducing the upper critical field linearly extrapolated to zero temperature. 



dT 

and the dimensionless reduced magnetic field, 



$0 



(15) 



B H 

we obtain a linear approximation for the critical line about T^, 

With the help of the dimensionless variables introduced so far it is easily worked 
out that the influence of the external magnetic field is reduced to discretization of 
the dimensionless in-plane kinetic energy, 

x^h{2n+l) (18) 

and the integrals of an arbitrary function / with respect to x are converted to sums, 

pc "0-1 

/ f{x)dx^2hy^ f{h{2n+l)). (19) 

In fact. Max Planck discovered the quantum statistics of the black-body radiation 
using the same replacement. Applying this procedure to the previously derived free 
energy at zero magnetic field, Eq. (12), we obtain 

F{e)^F{e,h) = -^, (20) 
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This expression represents the starting point for all further considerations. As a first 
step we address in the next section the Euler-MacLaurin method and its application 
to the sum over the Landau levels which appears in Eq. (21). 



2.2. Euler-MacLaurin summation for the free energy 

Near the critical temperature, when e <C c, one can consider formally c — > oo, and 
nc{h) « c/2h 00. Within such a local approximation the previous finite sums are 
transformed into infinite ones, 



2h /(e + /i(2n + 1)) = Sem / f{x)da 

n=0 { 



(22) 



where 



Of- 
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31 , 
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(23) 



is the Euler-MacLaurin operator for summation of series, in which we employ the 
Riemann and Hurwitz zcta functions, respectively. 



C(^) 



c(t/,i), c(i^,^) = E 



1 



n— ^ ' 



(24) 



instead of Bernoulli numbers. The Sem operator can be easily obtained exploiting 
the exponential representation of the standard translation operator T, whose action 
is defined as follows 



(25) 



If summed up from zero to infinity the above expression would give an infinite 
geometric progression, 

d 



1 



l-exp(6^)- 



n=0 n=0 L ^ 

Let us introduce now the fluctuational part of the heat capacity, 

1 92 

Using this physical observable and Eq. (21), one can extract the magnetic field 
dependent part of the free energy, 



(26) 



(27) 
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(28) 
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It is then straightforward to calculate the magnetization M and the nonlinear sus- 
ceptibility defined as x(e, h) = M/H, i.e. 

For the Meissner-Ochsenfeld (MO) state, for example, = — 1. Importantly, 

Eq. (29) incorporates the regularized free energy which, by virtue of Eq. (28), does 
not contain zero magnetic field part, 

j;eg(e, h) = F{e, h) - F{e) := R^g^uF{e) (30) 

where 



RegEM = Sem - i = 2/1 2 - / dx (31) 



Tl=0 

is the Eulcr-MacLaurin rcgularization operator. This method was applied"'^'^'"'^^'"'^^ 
for calculation of zero- field limit of the magnetic susceptibility, cf. Ref. 2. Hence, 
inserting the ft,^-term of Eq. (28) into Eq. (29) one finds^^ for 7/ ^ 

For illustration, let us analyze how this relation between susceptibility and heat 
capacity can be applied to the LD model. For arbitrary multilayered structure we 
can calculate the curvature of the lowest dimensionless energy band in c-direction 
^1(6*) = Sci{pz)lao, 

r^^^MO)\g=o- (33) 

According to this definition r parameterizes the effective mass in c-direction rric for 
an anisotropic GL model, \pz\ <C irh/s, 

It is now easily realized that the identity holds true, 

r={^)\ (35) 

which for A'' = 1 is the LD parameter r that determines the effective dimensionality 
of the superconductor, cf. the review by Varlamov et al}"^ In many other studies, 
e.g. Ref. 7, the wave vector k = p/?i has been used as well. In terms of the latter, 
for the dimensionless kinetic energy in the long-wavelength approximation we have 

[e.,{m.)+e.,{hk.)] ^ ^ ^ ^ ^3g^ 

do 
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where ea6(/ik) represents the in-plane part of the kinetic energy. Let us mention 
that the LD model is not only applicable to single layered cuprates with A'' = 1, 

but is it also to bi-layered cuprates {N = 2) in the limit cases 71 ~ 72 as well as 
in the case 71 <C 72 when formally = 1. That is why we use in our formulae an 
effective periodicity of the LD-model 

Seff = (37) 

For completeness we list below without deriving some of the well-known results 
within the LD-model. The single energy band has the form 

SoiP.) = aoo^i (g^) = ^^^^^ (1 - cos ^) , where (38) 

1 6 
ui{e) = -r{l-cose) = rsm'^-, (39) 

being parameterized by the Josephson coupling energy 



fi2 



"^^ = = m^s/N)^ ' r = 471 = {2NU0)/s f • (40) 
For the heat capacity one has 

and the magnetic susceptibility for a weak magnetic field applied in c-direction, 
according to Tsuzuki^^ and Yamayi,^° reads as 

^LD). X _ k^T , N 1 1 _ 1 Mo 1 

6 <pQ s y/e y/e + r b Hc2[v) [e + r) 



where 



$0 



(43) 



Before proceeding we feel it appealing to make some technical remarks concern- 
ing the representation of the general formulae for fluctuations in arbitrary layered 
superconductor. To be specific, we shall demonstrate how the expressions for the 
magnetic susceptibility, Eq. (42). and heat capacity, Eq. (41), within the LD model 
can be obtained as special cases of a general procedure described in the next sub- 
section. 



2.3. Layering operator L illustrated on the example of paraconductivity 

In the general formula for the density of the free energy, Eq. (12), the energies related 
to motion in c-direction ecj{pz) enter the final result solely via the fragment e+uij{6). 
Thus, in all such cases one can first solve the corresponding 2D problem and then 
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for a layered superconductor the result can be derived by merely averaging the 2D 
result with respect to the motion of the fluctuation Cooper pairs in perpendicular 

to the layers direction. Formally, this method reduces to introducing a layering 
operator L acting on functions of e; e.g. for the conductivity one would have the 
relation 

.(e) = La(-)(e) - ^ E / (e + ^f^ W) • (44) 

i=i 

In terms of the so introduced operator L the expression for the free energy, Eq. (21), 
takes the form 

F{e)=Foj\ij2 Lln(x + a;f )(^) + e) , Fo ^ ^^^7- (45) 

Besides for thermodynamic variables this operator works also for the fluctuation 
in-plane conductivity. Conforming with the work of Hikami and Larkin,^^ for the 
conductivity within the LD model we have to integrate the 2D conductivity with 
respect to the Josephson phase, 

a(e) = 1^"^° (e) = j ^u^^^^ (e + \r{l - cos 6)^ . (46) 

Given a system with independent 2D layers, having density in c-direction N/s, for 
zero magnetic field we have to average the well-known Aslamazov-Larkin expression 
for the static (zero-frequency) conductivity, 

where i?QHE = = 25.813 kfi. A simple integration gives 

. / X tCld)! f d0 1 1 

/.o(e;r).L - = f ^ J^T^^^^^^ - ^ 

We note that this integral determines both the heat capacity and magnetic sus- 
ceptibility for the LD model and is widely used for fitting to experimental data. 
Another important integral is 

= (2e + r) ln(\/I + \/f+7) - ^/e{e + r) - ln(4) [ , (49) 

which is used in representing the free energy at zero magnetic field, Eq. (12), cf. 
also Eqs. (129) and (130) below. Further, the e-derivative of this equation. 



1- = / ^ In + (1 - 00s 9)) = 21n( v^ + f^ ) 



(50) 



12 T. Mishonov & E. Penev 



is important for the calculation of the fluctuation part of the entropy and the 
density of fluctuation Cooper pairs. We provide also two other integrals employed 
in calculating the magnetoconductivity^^ 

.(LD)1_ WLD)gl_ e+^r 



dee [e(e + r)]3/2' 



. (LP) 1 ^ 1 I (LP) 1 ^ (e + r)£ + 

e3 2ae2 e + r)]5/2 " ^ 

Analogously, in the Maki- Thompson bi-layered model^° (A^ = 2) one summation 
should precede the integration/^ 

a{e) = L^^^V(2°)(e) = ^ ^ ^ [a^'^) (e + M^)) + '^^'''Ke + ^9))] , (53) 

i.e. in order to calculate the conductivity^^ and susceptibility^^ we have to add the 
terms 

1 , 1 2c + {^,+^,) ^^^^ 



In this expression both wi + u!2 and u)iu)2 are rational, cf. Eq. (4), and the integral 
(53) is reduced to the integral (48), 

- (MT) 1 e + 7i + 72 



e v/e [e + 2(7i + 72)] (e + 271) (e + 272) 



6 -|- —TW 

^ ^ fMT{e,h;r,w), (55) 



(e^ + rwe) (e^ + rwe + jr^w) 



where 



_ ^ 2 Jmax 7max 1 / „ , 7l , 72 \ ^ ( o i i "'^ 

r = 4^^ u=- = , w;=- 2H H = - 2 + u + - 

:^ + :^ Jmin 7min 4 V 72 7i / 4 V u 

Such a form, involving the w parameter, is convenient for fltting the experimental 
data since for both w = 1 and ^ 1 cases the LD approximation holds true, which 
is often found to give a satisfactory explanation of the experimental observations. 
For more detailed discussion the reader is referred to Ref. 21. The inverse relations 
for the above introduced parameters read as 

; 1 / 1\ r 1 r 

U={2w-l)+2^/w{w-l), Tmin = 2 ( 1 + - 1 41 7max = 2 (1 +'^) 41 (^6) 



and 



1 r 
2m„Tii;(0)2^'^"^4- 



Jmax = ao7max = ^ ^2 /r.\ o + 7" i^"^) 
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Importantly, at known effective mass niab the last equation gives the possibility one 
to determine the Josephson coupling energy between double Cu02 planes. 

Let us now illustrate in more details the action of the L operator. Towards 
this end we consider the famous Aslamazov-Larkin^^ formula for the 2D conduc- 
tivity cr^b°^(e), the results for the susceptibility —xfb'\^) due to V. Schmidt/^ 
A. Schmid^** and H. Schmidt, and Ferrell^^ and Thouless^'' fluctuation part of 
the heat capacity C^^'^^(e). With the help of the L operator for arbitrary layered 
superconductor, cf. Refs. 25 and 26, these three quantities can be generally written 
as 

1 ( 2toKT, \ N . 1 
= ^ \-^) - 'e' 

-Xabie) = -^l^o^eaM- L-, (59) 

where the "a5" -subscript in Eq. (58) indicates that the conductivity is in the ab- 
planes, while in Eq. (59) it indicates that the vanishing magnetic field is perpen- 
dicular to the same planes. It is immediately apparent that the common for all 
these expressions function Le~^ cancels when calculating the x/C, cr/C, and a/x 
quotients. In particular, the temperature independent ratio 

^0 = ^Ca6(0) TT = const 61 

provides the best method for probing the time constant tq parameterizing the life 
time of the fluctuation Cooper pairs with zero momentum, 

r(e) = ^, |Vp=o(i)|' oc expf-^rl • (62) 



e V T(e) 

The To constant participates in time-dependent GL (TDGL) theory; see for example 
the reviews by Cyrot,^ Skocpol and Tinkham,'' and the textbooks by Abrikosov^ 
and Tinkham.'* Within the weak coupling BCS theory in the case of negligible 
depairing mechanisms"^''*'^'^'''^^'^^'^'''^'^'^^'^^ tq satisfies the relation 

(BCS), rp (BCS) ^ 1 

Tq,^ ^b-'c _ TT (BCS) _ TT ft .^(BCS)/x_!0 ^ 1 

I -8' "16«' ^ ^'^-^""le^e' ^^^^ 

where Tq^^-* = 2tq^'~^^'' is the relaxation time constant for the order parameter be- 
ing two times larger.** At the present experimental accuracy this BCS value agrees 
well with the experimental data for the layered cuprates. Thus, the above observa- 
tion led us to propose the dimensionless ratio^^ 

= :;iScs) = 32-(^ = = -^UiO)^^ = const, (64) 
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to be used for more reliable experimental data processing; any deviation of frei from 
unity should be interpreted as a hint towards unconventional behavior and presence 
of depairing mechanisms. Notice also that the BCS value tt/8 = 0.393 in Eqs. (47), 
(63), and (64) is extremely robust, being originally derived for dirty 3D supercon- 
ductors, and the tqTc product remains the same^^'^'' for clean 2D superconductors 
and is not affected by the multilaminarity. The general formula for the fluctuation 
conductivity of a layered superconductor in perpendicular magnetic field can be also 
rewritten via the layering operator and relative life-time employing the 2D results 
by Redi,^^ and Abrahams, Prange and Stefen^^ (APS), cf. also Ref. 16, 

0-a6(e, /l) 



^re'TTI" l-/APs(e,^), 
Ion s 



where, ioi e + h > 0, 
/APs(e,/i) 



2hl 



-Ml 



V'(l 



2h 



(65) 



(66) 



is an universal dimensionless function of dimensionless reduced temperature e and 
dimensionless magnetic field h. The functions 



-lnr(z), V^i)(z)^-^(z) = C(2,z) (67) 
az dz 



are respectively the Euler gamma, digamma, and trigamma functions. This general 
formula is often utilized to process the experimental data for the paraconductiv- 
ity. We provide also several useful asymptotics of fxpsi'^ih) in different physical 
conditions, 

2 



/APs(e,/i) ~ < 



1 



e + h 



— In 2 
2h 



= 4 



T-T,2{Hy 

2 



h > lei 



(68) 



4 de^ 



For the LD model, for example, the {h <C e)-asymptotics gives 
Eq. (52), 



Trel 



mi~s 



V^(^T7) 2 [e{e + r] 



e(e + r) + |r 
1572" 



16,35,17 



3^2 



according to 



(69) 



Note that the classical Aslamazov-Larkin result, Eq. (47), is recovered for r = 0, 
/i = 0, and frei = 1- In the practical application to layered cuprates, however, we 
need to take into account the nonlocality effects. In -approximation to the GL 
theory we have to subtract the part of the corresponding cutoff area in the 2D 
momentum space. Thereby, the fluctuation conductivity is given by the difference 



Uab{e, h; c) = Caab{e, h) = aab{e, h) - Uab{c +e,h)^ C7„6(e, h) - aab{c, h), (70) 
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where a cutoff operator C is introduced, and the approximation is valid for e <C c. 
Similarly, for the magnetization we have the same "cutoff" expression which appears 
when calculating the truncated sums over the Landau levels, Yla" = ~ ' 
or integral with respect to the dimcnsionlcss in-plane kinetic energy. 



dx = 



dx 



dx. 



(71) 



As a rule the GL theory allows for ultraviolet (UV) regularization — every ex- 
pression can be easily regularized in the local (c — > oo)-approximation. Therefore 
the energy cutoff parameter c is not viewed as a tool for UV regularization, it is 
simply an important and immanent parameter of the GL theory, being of the order 
c ~ 1. The cutoff procedure has been essentially introduced from the beginning in 
the GL theory. *° Unfortunately, for many superconductors systematic studies for 
determination of the energy cutoff parameter are still missing. Here wo suggest only 
the simplest possible interpolation formula within the LD model for e <C c. 



crab[e,h) = C/LD(e,r «Trei — — 

8 -Kqhe s Ion s 



1 



1 



(72) 



which takes into account only the fist nonlocal correction. This simple expression 
fits very welP® the experimental data for YBa2Cu307_5 and Bi2Sr2CaCu208. 

As a last example of the action of the layering operator L we consider the 2D 
frequency-dependent paraconductivity at zero magnetic field. Taking the general 
expression for D-dimcnsional GL model"*^ and performing carefully the limit D — > 2 
(note, that there is an omitted term in the expression for the 2D conductivity in 
Ref 41) we get for the in-plane complex conductivity 



(73) 



or in expanded notation for singe layered superconductor. 



2Tok^Tc/h 



it/2 



^1 



e+r sin 



\e+r 



Soff^QHE Jo ""/^ e + rsin (/) 

where we have for the dimensionless real and imaginary conductivity, Ci (0) = 1 



(74) 



c;i{z) = — 



zarctan(2;) — - ln(l + 2^) 



arctan(2;) — z + z- ln(l -|- z"^) 



V T2(y) 



y 



—V / 

71" Jo 



dy, 
y"^ — z'^ 



dy. 



(75) 
(76) 



As usual in the above Kramers-Kronig integrations V indicates that the principal 
value of the integral is taken. For computer implementation of the L operator we 
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have to verify that, for r 1, L is simply equivalent to an incremental operator 
for the spatial dimensionality, 

a(°+i)(e)« ta(°)(e). (77) 

In many cases the GL results for integer dimensionality arc well-known and we 
can derive a generalization for a layered system. For both the MT and LD models 
the integration in Eq. (73) can be easily programmed, so we have a useful for- 
mula for fitting of the ultra high frequency measurements of cr*^(w). The original 
explicit expressions derived from retarded electromagnetic operator by Aslamazov 
and Varlamov**^ are too cumbersome to be used by experimentalists. Hence, one 
may realize that the GL theory is not some phenomenological alternative to the 
microscopic BCS theory (this scorn, dating back to the beginning of fifties, is still 
living even nowadays among students). The GL theory is a tool for applying the 
theory of superconductivity for the important for applications, let us say "hydro- 
dynamic", case of low frequencies and small wave- vectors. For e <C r the frequency 
dependent conductivity a*fj (uj), having dimension (O cm)^^, from Eq. (74) displays 
3D behavior, while in the opposite case of e ^ r it shows 2D character. For thin 
films of layered superconductors with thickness dmm we have to calculate the 2D 
conductivity cr*^^^^ = dfiimC, while for single layered films of conventional supercon- 
ductors, for example, we have to substitute in Eq. (74) formally Sgs = daim, and 
certainly r = 0. 

Having analyzed in detail the action of the L operator, we developed practically 

all technical tools necessary to proceed our investigation of the thermodynamics of 
Gaussian fluctuations and fluctuation magnetization. 

2.4. Power series for the magnetic moment within the LD model 

We will calculate in this subsection the nonlinear susceptibility by substituting first 
into the free energy, Eq. (28), the heat capacity, expressed via the susceptibility 
from Eq. (32). Then, the formula for the magnetization, Eq. (29), gives 

x(e,/.) = 6^(-l)"-i-£ h ___j^(2n)/.^(n-i)_^ (78) 

n=l ^ ' 

Taking the LD expression for the susceptibility at zero field Eq. (42), calculating 
the derivatives with respect to e by means of the relation 

d^J_ ^ (2m -1)!! 1 

Q^m ^ 2™e™ ^/e' ^ ^ 

and defining the relative susceptibility as 



(80) 
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we obtain 



V ,(e h-r)-12T^-'^Ul ' \ (2n + 2)!C(2n + 2) 

71=0 ^ ' 



X 



^ (2m-l)!!(4n-2m-l)!! 



^-J^ m!(2n - m)!(l + r/e)2«-'" 



Although these series is found to be a solution to the problem of calculating the 
fluctuational magnetization, 

Mich) = ix^^Ae^ _ xUc + e,h:r) \ 

^' ' $0Seff [6Ve(r + e) 6 ^(c + e)(c + r + e) J ' ^ ^ 

for the physical conditions of interest, i.e. an observable effect of magnetic field on 
the susceptibility, one needs to extend the series summation onto arguments /e^ 
beyond the radius of convergence. Analogous series has been already reported for 
the 3D paraconductivity. ^'^ One of the best devices for extending the convergence of 
series and also for calculating slowly convergent series is the e-algorithm*'*'*''' based 
on Pade approximants.'*'^ In the next section we describe a simplified version of this 
algorithm suitable for computer implementation. 



2.5. The epsilon algorithm 

The epsilon algorithm is a method for finding the limit L of infinite series 

N 

L = lim Sn = lim > a^, (83) 

i=0 

in case where only the first + 1 terms Ui, i = 0, . . . N, are known. The algorithm 
operates by employing two rows. The first one, called here auxiliary j4-row, is 
initially set to zero, i.e. 

a[°1=0, a[°1=0, a[°1=0,... ^[^1=0. (84) 

The second one is sequential iS-row loaded in zero-order approximation with the 

partial sums of the series 

Sq^ = ao, sf'^=ao + ai, ... sf^^ = ao + ai + a2 -\ h ajv- (85) 

The above assignments, as indicated by Eqs. (84) and (85), constitute the initial- 
ization phase of the £-algorithm. The essence of the latter consists of filling in the 
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SO called e-table 



q[0] 



A 

On 



[1] 





4°' 


4°' 


...\ 






^1 


^2 


...\ 












^0 


^1 


^2 








^3 






^0 


^1 


=■2 




5™ 


^2 


5™ 






^0 


=■1 


=■2 





[0/0] 


[1/0] 


[2/0] ... 


[1/1] 


[2/1] 


[3/1] ... 



(86) 



where according to the standard notations [j/k] = Pj{z)/Pk{z)\z=i designates a 
Padc approximant having power j in the nominator and, respectively, k in the 
denominator.^'' 

Starting from the and S^°^ rows every subsequent row is derived by applying 
the cross rule (known also as the missing identity of Frobenius) . To be specific, for 
calculation of the fcth A-row we have to solve the cross rule equation 



A' 



[k] 



olk-l] 



I ... NorthX 
West East 
South 



(South - North) (East - West) = 1. (87) 
Likewise, for calculating the fc-th S row we have to apply the same cross rule 



[fe] 



o[fe-l] 



i+l 



North 
West East 
\South ... 



South = North + (East - West) \ (88) 

Having applied the algorithm wc get in the S'-rows of the e-table, Eq. (86), a set of 
different Pade approximants to the limit L. The ith term of the fcth A-xaw can be 
easily obtained by 

4'=! = + (S%\^^ - Sf-^^y^ , for i = 0, 1, . . . , iV - 2fc + 1, (89) 

but for practical implementation of the algorithm, we can omit the index of the 
approximation and to use only one auxiliary row, updating it each time, 

Ai := Ai+i + {Si+i - Si)-^ , for i = 0, 1, . . . , iV - 2fc + 1. (90) 
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For the fc-th S'-row, the z-th term reads as 

5f 1 = ^l^-^J + (A^tii - ' for i = 0, 1, . . . , TV - 2fc, (91) 

and can be updated in the same manner as described for the ^-row, 

Si := Si+i + {Ai+i - Ai)-^ for i = 0, 1, . . . , TV - 2fc. (92) 

In order to find an estimate for the limit L of the infinite series, two different 
empirical criteria can be implemented. In the first one, the e-table is scanned for a 
minimal difference |/S'|^^^^ — iS]*^"^']. The limit L is then given by 



mm 

i.k 



Sf-'l. (93) 



This minimal difference gives also an estimate for the empirical error of the method. 
In the second criterion the £-table is scanned for the maximum of the East A-row 
element, cf Eqs. (88) and (91), 



max 

i.k 



A\k\ 



(94) 



The reciprocal of the maximum auxiliary value gives in this case the estimate for 
the empirical error of the method. It is the second criterion that we have used in the 
FORTRAN90 implementation of the e-algorithm given in Appendix A. Therein we 
have also made use of pseudo-inverse numbers in order to ensure provisions against 
division by zero in Eqs. (90) and (92), 

1 f 0, for z = 
11/2, tor z 7^ 

For an illustration, consider the first approximation. In the beginning we have 
for the first A-row according to Eq. (89) 

4'' = [(«o + ai)-(ao)]-' = ;^, 4'" = -^, ^l^-i = (96) 

The first S'-row then reads 

S^o^=ao + ai + — ^-r^, Si^' =ao + ai+a2 + — (97) 

1/02 -1/ai l/a3-l/a2 

and for the last element of the S'^l-row we have 

5'j^L2 = ao + ai + 02 H h aN-2 + aiv-i + (I/oat - 1/aN-i)- (98) 

The above approximation <S']^'_2 to limit L is nothing but the well-known 
Aitken's A^-method, which gives an exact result for the geometric progression 

5]^l = l + g + a^ + ••• + g^ 4^' = 5W = 5W = --- = 5j^l2 = Y^, (99) 
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for an arbitrary q^l. This fact can rationalize the success of the e-algorithm when 
applied to weak magnetic field series expansion of susceptibility. In the Euler- 
MacLaurin summation, Eqs. (22) and (26), we have a hidden geometric progression 
of translation operators. 

As a rule divergent series do not exist in physics; 99% of the divergent series 
born by real physical problems can be summed up by some combination of the 
Euler-MacLaurin method and the e-algorithm and the reason is lies in the analyt- 
ical dependence of the coefficients on the index. In the Gaussian spectroscopy of 
superconductors, for example, it is necessary series related to asymptotic expansion 
of Eulcr polygamma and Hurwitz zcta functions to be summed up, but the same 
methods could be applied to many other physical problems. The solution often can 
be derived by less efforts than required to verify that a series is divergent accord- 
ingly some strict mathematical criterion. Nowadays the mathematical education in 
the physical departments is conquered by scholastic mathematicians. Alas, none of 



the students of physics knows what really happens when we press the sin key of 
a calculator. On the other hand this is a commercial secret of the manufacturer. 
The physicists do not even lightly touch the brilliant achievements of mathematics 
indispensable not only for the theoretical physics but for experimentalist to fit their 
data as well. This is the motivation why we, following the spirit of the century 
of enlightenment, present in Appendix A a simple FORTRANQO program illustrat- 
ing the operation of the e-algorithm. Certainly, fysics is phun,^^ being in part art 
cosa mentale,'^'^ and every new software cannot be foolproof, but there are methods 
which must be taken into account in every complicated calculation. 

2.6. Power series for differential susceptibility 

Having calculated the relative dimensionless susceptibility by employing the s- 
algorithm we can recover the usual susceptibility from the dimensionless one, 

X(e,/i) = Xrci(e,Mx(e). (100) 

In order to take into account the effects of nonlocality the cutoff area in the mo- 
mentum space should subtracted out from the susceptibility 

(e, h) = Cx(e, h) = Xrei(e, h) x(e) - Xrei(c + e, h) x(c + e). (101) 

Then we can easily find the magnetization 

M{H,T) = x^{e,h)H. (102) 
The calculation of the differential susceptibility 

x''''He,h) = (^) , (103) 



T 



where H = Hc2{0)h, gives an alternative method to determine the magnetization. 
Next we define a dimensionless relative differential susceptibility 

iiie,h,r) = x^''''Ke,h)/x{e). (104) 
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For this variable, using Eq. (81), we have the series 



n=0 

2n 

X 



(2n + 2)!C(2n + 2) / h 



,2 



E 



22n+l / 22"+l 7r2n+2 I g2 

(2m- 1)!! (4n- 2m- 1)!! 



m! (2n — m)! (1 + r/e)^" 

m— ' ' 



which, just as done in deriving Eq. (81), can be summed up by means of the e- 
algorithm. For instance, in the local GL limit we have for the magnetization 

I^H I'll 

M= / x^''''\T,H)dH = x{e)H,2{0) / kie,h')dh' 
Jo Jo 

= x{e)H,2{0) Xreiie, h)h. (106) 

For the further analysis, however, it is more suitable to introduce a dimensionless 
magnetization 

rn^-^ = -^l.M, Mo.^^. (107) 



Then, using the relation 



Xrci(e,/i) = x / ~<^.h')dh' (108) 







h 

the result for the dimensionless fluctuation magnetization takes the form 



m(e,/i) = ^^=L=Xrei(e,fe)fe=^ r^—^ I k{e,h')dh' . (109) 
D -^e(e + r) o ^e{e + r) Jo 

In this section wc have calculated the magnetization by means of power series 
in the magnetic field assuming in the beginning h/e <^ 1. In the next section we 
develop another method for calculating the fluctuation magnetic moment which is 
appropriate for strong magnetic fields and allows for studying the high magnetic 
field asymptotics for large enough values of the reduced magnetic field, h/e^ 1. The 
overlap between these expansions about /i/e ~ 1 would be a test for the accuracy 
of the calculations. 

3. Strong magnetic fields 

3.1. General formula for the free energy 

In order to derive a general formula for the Gibbs free energy for arbitrary non- 
vanishing magnetic field we will start again by representing the free energy density 
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as a sum over the energy spectrum, Eq. (21), 

tlc-l 



F{e, h) = Fo 2h L ^ 



n=0 



In n + 



2h 



+ ln(2/i) 



where 



(110) 



(111) 



The first way to go in deriving convenient for programming formula is to calciilate 
the action of the L operator on the integrand, cf. Rcf. 48. In this case we write 
down the free energy as a finite sum over the Landau levels 



k T N 1 
h) = -jzklT^ — J2 + ^^(2/1 + 1)] , 

ra=0 



(112) 



where, according to Eq. (50), 

L^'^^^ ln[e + M2n + 1)] = ,,,V^ + H2n + l) + Ve + hi2n+l) + r^ ^^^3^ 

This formula is useful especially in the case of strong magnetic fields when the 
finite series are not too long. However, in order to have a good working expression, 
applicable to all cases, it is much better to solve the problem analytically. Towards 
this end consider the last term in the integrand of Eq. (110). The summation of 
this constant term and simply yields the cutoff parameter c 



2/i ^ 1 = {2h)nc = c. 

n=0 



Next we introduce a dimensionless function 



x{e, h) 



2h 



2h 



2H 



{T-T^H)) 



dT 



(114) 



(115) 



T=Tc-0 



which is the argument of some of the analytical functions we use in the following. 
Further, we have to present the sum in Eq. (110) as a difference of two appropriately 
regularized infinite series 



ln(n + a;) = Reg^ ^ ln(n + x) — Reg^ ^ ln(n + x). 

n—0 n—0 ric 

In fact, one does not have any other possibility except the C-regularization 



d 



-Reg^ X] + ^) = g^C(i^, z) 



n=0 



r(^) 



= In- 



(116) 



(117) 
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based on one relation between the Euler F-function and the Hurwitz i^-function, 
and the definition of the logarithmic function 



inz = lim = —z 

v—>o V av 



(118) 

i/=0 



According to the famous results by Riemann, the analytical continuations of the 
(■-function and the factorial n! are unique. Therefore the UV regularization of 
the partition function in the GL model in an external magnetic field is practically 
included in the second, Gauss definition of the F-function as a infinite product, see 
e.g. Ref. 49, 

/ t^-^(t-Ut = V(z)= lim -. (119) 

io ne^ooi;(2;+l)(^; + 2)...(^ + nc- 1) 

Let us recall some particular values, 

F(n + l)=n!, F(l) = 0! = 1, fQ^ = V^, (120) 

and the Stirling's approximation for Uc ^ 1, derived by Gaussian saddle point 
approximation applied to the first, Euler definition of the F-function, Eq. (119), 



-^j V^TTTic, In (ric!) « (ric + - j In nc - ric + In v27r. (121) 



For the local limit or for the case of weak magnetic fields we shall also make use of 
the asymptotic formulae for ^ 1 

\nT{z)^(z-^ynz-z+Un{2n) + ^^, (122) 



^(-)(.).(.-l)ln.-.+ ^, 



(123) 



2^-12^' ^'^(^)-^(2'^)-^2i + 6^- ^'''^ 

Substituting the Stirling asymptotics in the second Gauss definition, Eq. (119), and 
taking a logarithm we arrive at the function iIj^~^\z), generating the polygamma 
functions 

ij'^-'^\z) = Jhn^ |- ln(n + z) + (^Uc - ^ + In (ric) - ricj 



As a result, the above Gauss definition for lnr{z) solves the problem for UV reg- 
ularization of the infinite sum of logarithms, Eq. (117). The first derivative of this 
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equation gives the well-known definition of the digamma function ^p{z) = ip'^^^z), 



In particular, 

-V'(l) = CEuler 



n=0 



n=0 



n + 



lim 



^ - - In (nc) ^ = 0.577216 . 



(126) 



(127) 



All other polygamma functions are actually Hurwitz ^-functions with integer first 
argument > 2 and the sums are trivially convergent, 



i;{z) = i-i fm C(^ + l,z). 



(128) 



To summarize, wc have applied the well known ^-technique''''' for UV rcgularization 
of the partition function and revealed that the archetype of this powerful method 
comes from the century of enlightenment and finally we can bring the free energy, 
Eq. (110), to the form 



F{e,h) = -Fo |2ft L 
KT N 



inr 



1 e 



+ inr - + 



1 £ + C 
-I- 

2 



2h 



cln(2/i) 



4<6(0) s 



LC 



-(2/i)ln ^ 



27r 



;ln(2/i) 



(129) 



For weak magnetic field, /i «: e, cf. Eqs. (49), (50), (123), and (143), 



-{2h) In - e\n{2h) « -e [ln(e) - 1] + 

\/27r D e 



(130) 



This is our main analytical result and all thermodynamic properties now can be 
obtained via derivatives. However, having this analytical result it is trivially to check 
that it can be derived by finite sums. The latter do not require UV rcgularization 
and the Euler F-function is commonly available in many textbooks on mathematical 
analysis. 

3.2. Fluctuation part of thermodynamic variables 

Having an analytical result for the free energy we can easily find other thermody- 
namic variables by differentiating. The magnetization, for example, is given by the 
derivative 



M 



\dB 



1 



OF 



-Be2(0) \dh 



= -Mo m. 



(131) 
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where a dimensionless diamagnetic moment is introduced 



m(e, h) 



M _ c 
+ L 



inr(i±^Vinrr^±|±^ 

, , e + h + c 



2h 

2h 2h 



2h 



2h 



(132) 



In expanded notations within the LD model this formula, according to Eqs. (35), 
(37), (39), (43), and (46), reads as 



M(LD)(e,/i)^ 



- inr 



2h 



c + e + r sm 
2h 



Inr 



c + e + r sin 6 + h 



e + r sin (f) + h 
2h 



e + r sm 
2h~ 



c+e + r sin h 
2h 



e + r sm 



2h 



(133) 



where 



^9. For \e\, h r, c this general expression recovers the local 3D result, 
Eq. (180), analyzed later in Sec. 3.4, while in the opposite case of extremely high 
anisotropy r < |e|, /i <C c we get the local 2D result, Eq. (142). Here we want to 
emphasize the existence to mention a universal magnetization law at T = Tc, or 
e = 0, which can be observed for many high-Tc materials at strong magnetic fields 
/i > r 



_MK,B)^...i,...„(f^), 

where the universal function of the nonlocal magnetization 



(134) 



UM ^ ^ {inr(i + ^) - ^ + ^ [i - ^(^ + 0] } (135) 

is normalized so that Um{0) = 1, Um{oo) = 0, y = 2h/c. For conventional bulk su- 
perconductors the nonlocality effects on magnetization are well understood, see 
for example Refs. 51-56. To the best of our knowledge, the first observation 
of fluctuation-induced diamagnetism for a cuprate superconductor well inside the 
finite- magnetic-field regime was reported by Carretta et al.^"^ for YBa2Cu306+x- 
Soon after, analogous measurement was reported for Lai.9Sro.iCu04 by Carballeira 
et al.^^ Being familiar with the preliminary version of the present review (cf. Rcf. 59) 
Carballeira et al. have entirely based their interpretation and theoretical analysis on 
Eq. (133) and Eq. (142) below. Alas, we find it very disappointing and impolite that 
the authors of Ref. 58 do not give any credits (e.g. in the author list, acknowledg- 
ments, or references section) to the author (the first author of the present review, 
T. M.) of the theory they have used. We will not discuss in any detail their attitude 
and would instead refer to the Comment. 
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Returning now to the general expression for the magnetization, Eq. (132), we de- 
rive another expression for the relative differential susceptibility based on Eq. (109) 

dm 



K(e,/i)=6Ve(e + r) ( — 



= 6^/e{e + r) L 



+ 



4/i3 



2h^ 4/i3 ^ \ 2h 



€ + h + c 



2h 



(136) 



The comparison of this result with Eq. (105) is one of the best methods to check the 
accuracy of the programmed formulae. Analogously, differentiating the free energy 

with respect to the temperature T = (1 + e)Tc we derive the general formula for 
the most singular part of the entropy (neglecting the derivative of the T-prefactor 
in Eq. (129)), 



N 



4<fc(0) s 



e + h 
2h 



- t' 



2h 



(137) 



and the most singular part of the heat capacity 



1 d^F 
N 1 , 



477^2^(0) s 2h 



(1) 



2h 



(1) 



2h 



(138) 



This expression for C can be directly derived from the starting formulae (21) and 
(110). The sums for the heat capacity are convergent, cf. Ref. 40, and do not 
require any regularization. The simplest way to reproduce the analytical result for 
the free energy density, Eq. (129), is to integrate two times the result for its second 
derivative, i.e. that for the heat capacity, cf. Ref. 66. In general finite sums from 
to ric — 1 for logarithms and powers can be found in many textbooks on mathematics 
and all our results can thus be easily checked even by experimentalists. 

The fluctuation part of the entropy S is proportional to the mean square of the 
order parameter "t, i.e. the volume density of fluctuation Cooper pairs. The ther- 
mally averaged density deserves a special attention because it is the main ingredient 
of the self-consistent treatment of the interaction of order parameter fluctuations. 
This Hartree type approximation due to UUah and Dorsey^^ will be briefly described 
in the next subsection. 

In the following, for completeness, we will derive the local 2D asymptotics ap- 
plicable for \e\,h <C c. The substitution of the first term from Eq. (123) into the 
general formula for the free energy, Eq. (129), gives 



~ _ F{e,h) 



-(2/i)V 



e + h 
2h 



e\n{2h)+A{c)e + B{c) + 0{l/c), (139) 
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where for e -C c, cf. Eq. (12), 

rc+e 

fc{e) = A{c)e + B{c) Pi Inx dx ^ elnc + c{lnc - 1). (140) 
Jo 

This irrelevant for the fluctuation phenomena Hnear function of e gives constant 

additions to the free energy Fc = FoB{c), and entropy Sc = —FoA{c)/Tc and can 
be omitted hereafter. The subtraction of i^o/c from the free energy, Eq. (129), can 
be considered as a cutoff procedure for UV regularization. 



Reg|^i^(e, h) = F{e, h) - (F, - T,S,e), 



(141) 



which, when apphed, allows the analysis of the local GL approximation to be carried 
out simply as (c — > oo)-limit. Now a trivial differentiation gives for the dimensionless 
magnetization, being a positive quantity. 



m2D(e, h) 



-M{e,h) _ _e_ 
" 2h 



Mo 



1 5/2D 
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2h 



1 



(-1) 



2h 



(142) 



2 dh 

This result is also a local (c — > oo)-asymptotic of Eq. (132), which for /i <C e yields 

m2D » h/%e. (143) 
In the general case the local approximation gives m = Lm2D, or for the LD model 



m(e, h;r) = — 



M 
Wo 



Jo ^/2 



■ r sm 



2h 

-inr 



r sm 



2h 



+ ■ 



1 



e + rsin (j) 1 
2h ^ 2 



- 1 



ln(27r) } . (144) 



The next differentiation with respect to the magnetic field, using Eqs. (42), (103), 
gives the relative dimensionless siisceptibility 



ee(§).:.(^)n.-4,<.,(^,^ 



(145) 



which is also a local c ^ h, \e\ asymptotic of Eq. (136). For the LD model after 
averaging with respect to the Josephson phase, according to Eq. (39), we obtain 



K{e,h;r) = 12 



r/2 



it/2 



e + r sin 

2/r~ 



1- 



e + r sin 

2/r~ 



C 2, 



e + r sin 
2h~ 



(146) 



This final result can be directly compared to low field series expansion Eq. (105). 
Similar differentiations of the free energy, Eq. (139), with respect to the tempera- 
tures gives the most singular part of the entropy 

'e + /i' 



S2D 



d ~ 



2h 



ln(2/i) 



Tc5(e,/i)/i^o, 



(147) 
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and of the heat capacity 



LC2D = -L^/2D(e,/l) = ^V(^) 



e + h 
2h 



= T,C{e,h)/Fo. 



(148) 



Restoring the T prefactor instead of in Eq. (12), as was done in Eq. (129), we 
arrive at a shghtly different expression for the fluctuation part of the free energy 
F = Fq{1 + e)/2D(ej h) and the heat capacity 



tc2D =-(l + e)^(l + e)t/2D(e,/i) 
= (1 + e) [(1 + e) Lc2D + 2 Ls2D 



which gives 



C(e, h) = 



4<fc(0) s 



1 + e 
2h 
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€ + h 

2h 
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,(0) 



e + h 
2h 



21n(2/i) 



(149) 
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For zero magnetic field we have 

(1 + e)K N 



C{e,h = 0) = 



4<\(0) s 



(1 + e) L ^ - 2e L In ^ 



which in the LD model takes the form 



+ 



(151) 

-2e21n^ — . (152) 

These expression differs from Eqs. (41), and (60). However, the (1 + e)^ w 1 + 2e 
correction and the less singular part of the heat capacity 2(1 + e)Fo Ls2d/Tc, which 
appears due to differentiation of T in the numerator of Eq. (12) and Eq. (129), are 
difficult to be identified experimentally. 

For the superconducting phase below the critical temperature, < — e ^ 1, one 
has to take into account more or less space homogeneous order parameter ^/^ which 
minimizes the nongradient part of the free energy density F = a(e)ng + , 



*e = \/ao{-e)/b, Tie = = ao{-e)/b. 
The fluctuations around this minimum 



//^2 



* = *e + + n = ^^=n, + 2*,*' + (*') + (*") 



(153) 



(154) 



should be considered as a small perturbation, thus only the quadratic term in the 
free energy is taken into account. 



F{c < 0) = a(e) + ^bn' « -^a^e^ + ao{-2e) [l (M/')' + (M/")' 



(155) 
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The first term in this equation corresponds to the jump in the heat capacity AC = 
a^/bTc at T^. The linear term cx ^' simply cancels. The phase fluctuations a 
('I'") arc coupled to the plasmons and vortexes but they are irrelevant for the 
thermodynamic fluctuations signiflcantly below T^. In this way mainly fluctuations 
related to the modulus of the order parameter are essential for the heat capacity 
below Tc- Finally, the comparison of the second term oc ^' in Eq. (155) with the 
corresponding expression above Tc 



F{e > 0) = a(e) + -bn^ w ao(e) 



(156) 



provides a prescription to derive the fluctuation part below Tc from the fluctuation 
expression for the normal phase above 
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2 (-2e) 



l1. 
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(157) 



Applying this prescription to Eq. (152) results in the following expression 



^ ' ^ s 2 



+ 



- 2e 2 In 



v'(-2e)(r-2e) ./(Z^ + 



(158) 



These fluctuation part as well as the the phonon heat capacity should be subtracted 

from the experimental data in order to extract the jump AC and related to it 
penetration depth Xab{0)- Such a procedure, in fact, gives a purely thermodynamic 
method to determine the latter quantity. 

The dimensionless functions Eqs. (139), (142), (145), and (148)) derived with the 
local approximation are just as important for the thermodynamics of the layered 
superconductors as is the APS function for the paraconductivity, Eq. (66). The 
operator L gives the possibility to extend the 2D analytical result for layered or 
even isotropic 3D superconductor. Additionally the C operator gives the energy 
cutoff approximation for the nonlocality effects in the conducting Cu02 planes. 
Therefore the analytical 2D result plays a key role for the fluctuation phenomena 
in layered superconductors. 

We will finish the analysis of the local c » \e\,h 2D approximation r, i.e. 



2um\ 

s J 27r^2^(0)' 



(159) 



with the important case of strong magnetic field h |e|. Under these conditions 
(|e|, r ^ /i -C c) the layered superconductors display a magnetization corresponding 
to the local 2D one in strong magnetic fields. The substitution of e = in Eq. (142), 
using Eq. (120), recovers the result by Klemm, Beasley, and Luther^^ 



7h{h » r, e < /i) fa 0.3465735902799726 . 



-M 



In 2 KTcN 



(160) 
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In the concise review by Koshelev^^ on the properties of 2D GL model the cal- 
culation of 5 In 2 Ri 0.346 by infinite series with three decimal digits accuracy is 
described in great details. 

3.3. Self- consistent approximation for the LD model 

The bulk (3D) density of the fluctuation Cooper pairs n(e, h) can be calculated 
from the general expression for the Gibbs free energy Eqs. (8) and (12). The dif- 
ferentiation with respect of the "chemical potential" of Cooper pairs ^cp = — cioe, 
according to the relation A^cp ~ {dG/diJ,cp)T,H, gives 

n{e,h) = -( \^nf) = -^F{e,h) = -^S{e,h;c). (161) 

This formula can be alternatively derived by summation of the Rayleigh-Jeans 
asymptotics of the energy distribution of the fluctuation Cooper pairs 

^' ' V ^ .ej{p,p^) + a s J (27rfi)2p2/2ma6 + aoe' ^ ' 

P'P"^ IpKpc 

see for example the monograph by Patashinskii and Pokrovsky.^^ 

Let us give an illustration for zero magnetic field. In this case for the density of 
fluctuation Cooper pairs, using Eq. (Ill) and Eq. (50), we obtain 

n{e, 0) = — 2 In . (163) 

This formula sets the stage for the self-consistent treatment of the order parameter 
fluctuations in the LD model in which the nonlinear term is replaced by its average. 
The idea has its origin in the Maxwell consideration of the ring of Saturn; probably 
it is the first work on collective phenomena in physics. Having no possibility to 
consider motion of all particles in detail we must search for some approximation. 
Within a self-consistent picture, the motion of every particle creates an average 
potential in which the others are moving. From the dust of the ring of Saturn 
to the Cooper pairs in cuprates the idea is the same, only the mechanics slightly 
changes. In the self-consistent approximation the nonlinear term in GL equations 
gives an addendum to the linear one 

arcn(e, /t) = «oe + & « ( ) , (164) 

where the coefficient h = hN j s can be expressed via the jump of the heat ca- 
pacity AC at the phase transition or, which is more convenient for the high-Tc 
cuprates, via the extrapolated to zero temperature penetration depth 1/A^(,(T) = 
Mon(T)e*V"^a6, n(T) = -a{T)/h, 

^=^ = ^^'{^^) ' ^^^^=8^(a„,(o!L(0)) ' ('^^^ 



Thermodynamics of Gaussian fluctuations and paraconductivity in layered superconductors 31 



where k^^ = Xab{^) / ^ab{0) is the GL parameter. One can easily check that Eq. (165) 
has the same form in Gaussian units, where fj,^''^^^^^ = Att. Introducing the renormal- 
izcd reduced temperature Cren > for the normal phase we have the self- consistent 
equation 



. T b . 

Cren = m — H n(eren, h), 



(166) 



where n(e, h) is calculated by means of Gaussian saddle point approximation.^^ For 
the LD model this equation, by virtue of Eq. (163), takes the form 



, T + Eren + \/c + Cren + r , -> (LD) C + Eren 

Eren = lu — + e^,, 2 In . — == = e + e^j L In , 

J-c Y Cren ~r V^ren "r 7* Cren 

where the dimensionless parameter 



(167) 



$0 



4^,(0) s AC 



(168) 



is closely related to the Ginzburg number; cf. Eq. (60) which now reads 

and the review article by Varlamov et al}"^ At T^, for e^, <C r <C c, Eq. (167) gives 



4c 



(170) 



and the effective heating AT = Tceron.c constrains the fluctuation variables at T^.. 
To provide an order estimate we take for illustration Seff = 1 nui) Aab(O) = 207 nm, 
Tc = 100 K, ^ah(O) = 2.07 nm, k^l = 100, = 1.381 x IQ-^^ J/K. The substitution 
of these values in Eq. (168) gives 



87r^ X 1.381 
1000 



11%, 



6^2 



2 X 10" 



(171) 



In the case of nonzero magnetic field the self-consistent equation for the rcnormal- 
ized reduced temperature, Eq. (167), according to Eqs. (137), (161), and (166), 
takes the form 



T 

In— -Fee., L 



2h 



2h 



or, within the LD model, 



(172) 



'eren + h + r sm 
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(173) 
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cf. also Ref. 65. For weak magnetic fields, /i -C e, using the asymptotic formula for 
the digamma function, Eq. (124), we recover Eq. (167). The formulae pointed out 

could be easily programmed for the sclf-consistcnt LD fit to the paraconductivity 
near to the critical temperature T^. With the foregoing discussion we finish the 
analysis of the thermodynamics of layered superconductors. We only note that all 
final formulae can be used to fit the experimental data. Before proceeding however, 
for reliability sake, it is necessary to check if the formulae implementation correctly 
reproduces the 3D limit case r ^ oo. 

3.4. 3D test example 

Every layered superconductor near the critical point |e|, h r displays 3D behavior. 
For high-Tc cupratcs, however, r <C 1 and 3D behavior can be observed only in 
crystals of extremely high quality. Due to fluctuation of the stoichiometry and 
of the Tc 3D regime of Gaussian fluctuations may not occur. However there are 
many conventional layered compounds with moderate anisotropy, r < 1, to which 
the 3D behavior has broader applicability. The 3D case can be derived as (r 
oo)-asymptotics if the parabolic band approximation loi{6) « r^^/4, Eq. (34), is 
substituted into the L operator, Eq. (46). Using the variable x(e, h) from Eq. (115) 
and a new dimensionless variable q, defined as 

1=\H'^^ de = 2Mdq, (174) 

y 8h h 8h V r 

we get for the regularized sum of logarithms in Eq. (112) the local approximation 

t(LD)-r- / 1 e\ ^ f2h ^ f+°° dq ^ , 2^ 

^Reg,^ln n + - + - «2^/-Reg,^ / ^ln{n + x + q^). 

n=0 ^ ' 71=0"'-°° 

(175) 

The UV regularization in this expression is carried out with the help of the equation 
C E /7 ^ + ^ + 9^) = C (- ^, , (176) 
which can be easily proved using derivatives of ^-functions 

^C{u,x) = -uC{iy + l,x). (177) 

The second derivative of Eq. (176) is trivially convergent; the essence of the (- 
function regularization lies in the omission of an arbitrary linear function A{c)x + 
B{c), being analytical with respect to e and therefore irrelevant to the critical 
behavior, cf. Eq. (141). In fact c ~ 1 but having dropped A{c) and B{c) we can 
obtain the local approximation, |e|,/i <C c, as c — > oo even if A{oo) = oo and 
B{oo) = oo. The substitution of this UV regularization in Eq. (112), using Eq. (35), 
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gives the result by Mishonov^^ for the fluctuation part of the Gibbs free energy 



This result was confirmed by Baraduc et al.,^^ using the same notations, with the 
C-function presented implicitly. 

In order to bridge the 3D rcsidt with the notations introduced for layered systems 
we can rewrite the coefficient in Eq. (178) as 



^ 

27r ^a(0)6(0)€c(0) 



(179) 



Now differentiation i^(e, K) with respect to the magnetic field we obtain for the di- 
mensionless magnetization, in agreement with the result by Kurkijarvi, Ambegaokar 
and Eilenberger^^ 



m(e, K) 



31? 
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2' 2 "^2/1 



3'' U' 2 ^ 2h) 2h 



(180) 



The subsequent differentiation with respect to the magnetic field gives the difi^eren- 
tial susceptibility. In the particular case of strong magnetic fields, e -C /i, the local 
approximation to the GL model, Eq. (180), gives the well known-result by Prange^^ 
with an anisotropy correction multiplier^^ ^ab{0)/^c{0) 

m(0, h) = 3V2 x 0.0608885 ^, -M = ^ V^, 

(181) 



where for the values of the C-function we have 
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2' 2 
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■1 ' r- 
\/2, 

Zeta[-l/2] 



-ij = Zeta[-l/2,l/2] = 0.0608885..., 
-0.207886 .... (182) 



The syntax Zeta[ ... ] is used in the commercial software Mathematica® .^^ We 
stress, however, that these are only test mathematical asymptotics for c ^ 00. 
For the magnetization, as well as for every quantity exhibiting UV divergences in 
the local limit, the nonlocal effects are strongly pronounced simply because the 
contribution of high momenta is significant. That is why the local approximation 
could be quantitatively fairly good for fitting to the data for fluctuation conductivity 
and heat capacity. For the magnetization in strong magnetic field regime wc have 
to take into account the effect of nonlocality by fitting the energy cutoff parameter 
. A systematic procedure for determination of the parameters of the GL theory 
is developed in the next section. 
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4. Some remarks on the fitting of the GL parameters 



4.1. Determination of the cutoff energy 

Let we start with designing a general procedure to fit some parameters of the GL 
theory which employs only data for the in-plane paraconductivity. Later on we 
shall address the advantage of investigating several variables simultaneously. The 
first step is to extract the fluctuation part of the conductivity from the temper- 
ature dependence of the resistivity R{T). For layered cuprates the resistivity of 
the normal phase is to within good accuracy a linear function of temperature, 
Rn{T) = An + BbT, and we can fit the coefficients An and Bn far enough from 
the critical temperature T^, e.g. in the temperature interval (1.5 T^, 3 T^). After that 
we can determine the experimental data for the fiuctuation conductivity 

ai = R{Ti)-^ -{AR + BnTi)-^ (183) 

for all experimental points i = 1, . . . , A^exp- For bi-layered cuprates, such as YBa2Cu307_5 
and Bi2Sr2CaCu208, one can attempt fitting the data with the formula for the bi- 
layered model, Eq. (55), where an arbitrary life-time frei and cutoff parameter c are 
included in the interpolation 

e2 N 

fT(e; frei, r, w, c) = J^'^reiy [/MT(e; r, w) - fuTic + e; r, w)] , (184) 

where 

f , ^ e+^rw ^(MT) ^(MT)1 

jMT{e;r,w) = == L /Aps(e, /i = 0) = L 

J (e^ + rwe) (e^ + rwe + jr^w) ^ 

(185) 

As a next step, if necessary, one may fit the data using logarithmic plot that gen- 
erates the dimensionless deviations from the ln(tTa6(T))-values, 

Xi{r, w, c) = In [/MT(ei; r, w) - /mt(c + ej; r, w)] - In I | . (186) 

For the Xi data we can calculate the mean value, the averaged square 

^ ^ cxp ^ ^ exp 

and the dispersion 

S{r,w,c) = {x'^) -{xf. (188) 

The fitting procedure is then reduced to numerically finding the minimum of the 
dispersion 



S{ro,WQ,co) < S{r,w,c) 



(189) 
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in the space of parameters {r,w,c). We have to start from some acceptable set 
of parameters, for example, tq = wq = 1, and cq = 5, and to search for the 
minimal value in certain range, e.g. ro =G (0, 1), wq E (1.30) and co = ^ € 
(0.2,2). It is possible that the parameters of the normal resistivity be corrected 
by the same procedure for minimization of the dispersion S{r,iu,c,Aji,Bf{). The 
contemporary methods of the mathematical statistics, such as the bootstrap and 
Jack-knife, can tell us how reliable is the set of the fitted parameters; the simplest 
possible realization is to decrement sequentially iVexp by one and to investigate 
the distribution of the fitted GL parameters at every step. For example, the w 
parameter is almost inaccessible since for w = 1 and w ^ 00 wc have LD-type 
temperature dependence of the paraconductivity. On the other hand if we try to 
fit the paraconductivity far from the critical temperature, e.g. T G (1.02 Tc, 1.15 Tc) 
we can easily find some estimate for the cutoff parameter c. In any case a good fit 
would be useful because as a by-product we determine the life-time of the fluctuation 
Cooper pairs 



frei = exp {-(x}{ro, Wo, Co)) . 



(190) 



The same procedure can be applied to the magnetic susceptibility at vanishing 
magnetic field which, according to Eq. (61), is proportional to the conductivity, 
or for the susceptibility in the LD model which, according to Eqs. (42), (72), and 
(168), reads 



where 
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+ r) V(c + e)(c + e + r) 
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-5- - 27r/xo-p-U(0)— - 



The general formula for the conductivity, Eqs. (58), (70), 

o-ob(e,«) = g-^^^ L C/APs(e,ft), 
which for single layered superconductor reads, cf. Eq. (133), 
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(191) 



(192) 



(193) 



, , 1 e + r sin 

^i2 + ^r- 



(194) 
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gives another possibility for determining the energy cutoff parameter c. As most 
appropriate regime we recommend that the measurements of the conductivity as a 
function of the magnetic field to be carried out at the critical temperature T = Te- 
la this case for strong magnetic field, h ^ r, the layered superconductors with 
strong anisotropy r <C 1 show 2D behavior. The substitution e = in Eq. (194) 
gives another universal law derived within the GL theory with energy cutoff 



1 B aab{e = 0,h) 

2 Be2(0) f,el(eVl6/lSeff) 2 '^^P'^' 



e To 



B 



C5e2(0) 



where y = 2h/c, cf. Eq. (135), 
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(195) 



(196) 



[/g.(0) = 1, and Ua-{oo) = 0. At best, the universal dimcnsionlcss conductivity 
Ua- oc Ba{B) and magnetization Um oc M have to be fitted simultaneously using 
the data for the same crystal and common dimensionless argument oc B. Similar 
universal scaling law for the heat capacity can be derived from Eq. (138) 



B C{T,,B) (2 B 

— Uc 



C(2,i)Bc2(0) e,.AC 



cBe2(0) 



(197) 



where 



Uc{y) 



C(2,i) 



(198) 



but the accuracy of thermal measurements is probably not high enough in order for 
this to be experimentally confirmed. 



4.2. Determination of the coherence length ^a6(0) 

The fit of every fluctuation variable as a function of the dimensional magnetic fleld 
h, the conductivity 

C7(e,/i) = CT(e) + ACT(e,/i), (199) 

for example, provides a method for determination of i3c2(0) and ^afe(O). At weak 
magnetic fields, /i <C e, the magnetoconductivity is proportional to the square of 
the magnetic field Acr(e, h) = cr(e, h) — a{e) oc B'^. For this small negative quantity, 
< -Aa{e, h) «: (j{e), the APS result, Eq. (68), reads^^ 

-Aa(e,ft)«^^a(e), (200) 
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where h = 27r^^^(0)i?2/$o- The common multipher tq from Eq. (65) is obviously 
canceled in this relation because, roughly speaking, the transport takes time even in 
the presence of magnetic field. Wc note that a multiplier fid was misintroduccd by 
M. V. Ramallo in Ref. 36 in the right-hand-side of the above equation (see Eq. (4) 
in Ref. 36). Thereby the old experimental data in Ref. 36 have been apparently 
processed by employing erroneous expression and therefore the discussion related 
to Fig. 2 in Ref. 36 is physically unsound. As a consequence, the life-time constant 
of metastable Cooper pairs in cuprates is still waiting for its first experimental 
determination. Nevertheless the novel theoretical result that the life-time constant 
To and the diffusion coefficient of the fluctuation Cooper pairs ^a(,(0)/To can be 
determined from the a/x-ratio remains unchanged. 

Returning to Eq. (200) we note that after two-fold integration of the relation 
(200) in some temperature interval, e.g. (ca, eb) = (0.03,0.09), the "noise" in the 
experimental data is already irrelevant and we can rewrite Eq. (200) as^^ 



U{0) = Ib 



J'' de' £ (-Aa(e", h)) de" 



o-(ea) - o-(e6) + (£6 - 



1/4 



(201) 



where Ib is the magnetic length 




25.6 nm 

~7Wy 



(202) 



For practical application we have to take into account that far from the critical 
temperature, even for T — = 15% the fluctuation conductivity is negligible 
(t(0.15) k, (t(oo) = 0. That is why in acceptable approximation we can take £{, = 
0.15 « "oo". For e > et the temperature dependence of the magnetoconductivity in 
the numerator of Eq. (201) can be an extrapolated LD fit. 

However, due to the strong critical behavior — Act oc h'^/e^ for e^ r the influence 
of the interval {ei,,oo) can be neglected. In such a way, after a partial integration, 
we arrive at a simpler equation for determination of the in-plane coherence length, 
cf. Ref. 36, 



1 



e' {-Aa{e', h)) de' - e / (-Aa(e", h)) de" 



1/4 



const, 
(203) 



where the integrations should be performed in the whole experimentally accessible 
temperature range above (1 -|- e)Tc. This result of the Gaussian fluctuation theory 
does not depend upon the ro parameter, effective mass of Cooper pairs ruab, and the 
space dimensionality. We consider this procedure for determination of the coherence 
length ^ab(O) as being the best one, as it is model- free and does not depend on the 
multilaminarity of the superconductor, i.e. on the dispersion of Cooper pairs in c- 
direction ec,j{Pz)- Equation (203) has the same form for both strongly anisotropic 
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high-Tc cuprates and bulk conventional dirty alloys. Of course, methods particularly 
based on the proximity to the critical line Hc2{T) can be very useful in determining 
Coh(O) especially in the case of strong magnetic fields. For example, Eq. (68) gives 
another appropriate formula 

TV 4 

aabie, h) « frei-777: (204) 

applicable for <C e + /i -C /i. Similar result, cf. also Eq. (142), 

k T N h 
M = -Mom^-%^-^= (205) 

$0 s ^(e + h)ie + h + r) 

can be derived under the same physical conditions from the formula for the fluctu- 
ation magnetic moment, Eq. (133), using the approximations for < a; <C 1, 

r(a;) K-lnx-^ ln(27r), i){x) w --. (206) 

The experimental investigation of the conductivity, Eq. (204), and magnetization, 
Eq. (205), is probably the best way to extract the upper critical field Hc2{T) for 
high-Tc cuprates; the Ha/M quotient near the critical line is 2/3 of the cr/x quotient 
for weak magnetic fields. 

4.3. Determination of the Cooper pair life-time constant tq 

Having a reliable estimate for the coherence length, the life-time constant of the 
metastable Cooper pairs above Tc can be determined via the a/x-quoticnt, Eq. (61). 
We believe that this method will become a standard procedure in the physics of 
high-Tc materials. Certainly the most transparent method is just the fit to the 
phase angle of high-frequency complex fluctuation conductivity 

(Pcr[coT[eren)) = arctau = arctan^; . (207) 

The state-of-the-art electronics gives such a possibility, but unfortunately the first 
experiments of the type'*^'^® was not performed in the Gaussian region. For the 
development of Gaussian spectroscopy which will give results relevant for the micro- 
scopic mechanisms of superconductivity we recommend the use of the conventional 
thin films and high-quality low temperature cuprate films, such as Bi2Sr2CaCu208. 

4.4. Determination of the Ginzburg number and penetration depth Xab{0) 

The applicability of the self-consistent approximation in the theory of fluctuation 
phenomena in superconductors is strongly limited by the quality of the samples. The 
fluctuation of the critical temperature ATc, e.g. due to the oxygen stoichiometry 
in cuprates, should be small enough, ATc ^a'^c, and this has to be verifled 
empirically. If the a/x ratio remains temperature independent for e < 3% and 
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both a{e) and x(e) demonstrate weak deviation from the LD fit obtained from 
the range e e (3%, 9%), this could be considered as a hint in favor of the self- 
consistent approximation. In this case e^,j can be fitted by substituting the solution 
ere„(e) of Eq. (167) into the LD fit to a'^^^^{e,en) and x^^^K^ven)- We note that the 
reliability in fitting is determined by the condition whether the self-consistent 
approach and the use of Cion significantly improve the accuracy of the fit to the 
experimental data near Tc- According to Eq. (168) we can parameterize e^,^ with 
the help of the penetration depth Aafc(O). In any case, an evaluation of such type 
should be a part of the complete set of GL parameters of the superconductor. 
Another possibility for the thermodynamic determination of the penetration depth 
Aa6(0) is provided by the jump in the specific heat at the critical temperature, 
Eq. (165). As a rule the accuracy of the determination of the penetration depth by 
the thermodynamic methods cannot be high, especially for liigh-Tc cuprates where 
the phonon part strongly dominates. An acceptable value of InK^^ derived from 
the heat capacity is necessary for the establishment of a coherent understanding of 
the superconductivity; there is no doubt that the direct investigation of the vortex 
phase of the superconductors or vortcx-frcc high-frequency measurements constitute 
the best methods for determination of Aab(O). 

5. Discussion and conclusions 

In the attempts to systematize the available results we had to derive in parallel new 
ones too. We shall summarize the most important of them starting with remarks 
concerning the theory. As the ultimate result we consider the representation of 
the fluctuation part of the Gibbs free energy by the Euler F-function Eq. (129) in 
Gaussian approximation. This result trivializes the derivation of all thermodynamic 
variables, such as fluctuation magnetization Eq. (133), or fluctuation heat capacity 
Eq. (150). To our knowledge this is a novel result, but we find it strange that it 
remained unobserved given the great attention which the fluctuations in high-Tc 
superconductors have attracted. The importance of fluctuations was mentioned 
even in the classical work by Bednorz and Miiller. Fluctuations in superconductors 
were among the main topics in many scientific activities; the F-function is well- 
known to all physicist; the mathematical physics behind the 2D statistical mechanics 
is well developed; polygamma functions can be found in a number of BCS papers, 
and finally the solution turns out to be on a textbook level. Just the same is the 
situation for the 3D GL model. In this case the solution for the free energy is given 
in terms of the Hurwitz i^-functions. Analogous result gave the name of one of 
the most powerful methods in the field theory — ^-function method for ultraviolet 
regularization, but this method was never applied to the most simple problem of a 
3D GL model related to numerous experiments in the physics of superconductivity. 

Another simple but useful detail is the layering operator L, Eq. (44), which 
allows us to extend the 2D result onto layered superconductors and even to 3D 
superconductors. The method can be applied not only to the thermodynamic vari- 
ables but to the fluctuation part of the kinetic coefficients as well. In this way we 
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obtained useful formulae for the in-plane fluctuation conductivity in perpendicular 
magnetic field, Eq. (194), and for the high-frequency Aslamazov-Larkin conduc- 
tivity in layered superconductors, Eq. (74). Wc proposed further convenient r-w 
parameters, Eq. (55), for the bi-layered model which could be utilized for experi- 
mental data processing of the fluctuation phenomena in bi-layered cuprates, such 
as YBa2Cu307_5 for example. 

The representation of the thermodynamic variables via polygamma functions is 
very helpful at strong magnetic fields but due to the presence of the magnetic field 
in denominator these results cannot be directly applied to zero-magnetic-field limit. 
For small magnetic fields, on the other hand, we have to use asymptotic formulae 
for polygamma and C-functions with large arguments. This is the reason why the 
weak-magnetic-field expansion of the magnetization and the other thermodynamic 
variables has so bad convergence. In order to fit the experimental data for the 
magnetization in weak magnetic field using the new analytical result for the LD 
model, Eq. (82), we arrive at the problem for summation of divergent asymptotic 
series. At least for experimentalists this is a nontrivial problem which led us to give 
a prescription for usage of series from the theoretical papers. There is no doubt that 
the £-method is one of the brilliant achievements of the applied mathematics of XX 
century. However, it turns out that this method was not cast in an suitable form to 
be employed by users like experimentalists having no time to understand how the 
underlying mathematics can be derived. That is why wc presented an oversimplified 
version of this algorithm illustrated by a simple FORTRAN90 program. The latter 
can be also used for calculation of the differential nonlinear susceptibility at finite 
magnetic field, Eq. (105), which is another novel result in the present work. 

Let us now address the simple final formulae that can be directly used for exper- 
imental data processing. First of all we advocate that the relation between fluctu- 
ation conductivity and magnetoconductivity, Eq. (203), provides the best method 
(shortly announced in Ref. 36) for determination of the in-plane coherence length 
^afe(O) in layered high-Tc cuprates and conventional superconductor superlattices 
and thin films. Having such a reliable method for determination of ^a6(0), the 
Cooper pair life-time spectroscopy can be created^^ on the basis of determination 
of the life-time constant tq by the a/x quotient, Eq. (61). 

Usually science starts with some simplicity, thus it is surprising that the tem- 
perature independence of the a/x, xl^i ^md cr/C quotients has not attracted any 
attention in physics. The question of whether the high-Tc cuprates are BCS super- 
conductors, or they have a non BCS behavior, consumed more paper and brought 
more information pollution than that about the sense of life, about the smile of 
Mona Lisa. Now we possess a perfect tool to check whether this sacramental ^ 
BCS ratio, Eq. (63), still exists in the physics of high-Tc superconductivity. A care- 
ful study of the relative life-time constant tq by the a jx ratio, Eq. (64), will provide 
a unique information on the presence of depairing impurities in the superconducting 
cuprates. The doping dependence of this ratio will give important information for 
the limits of applicability of the self- consistent BCS approximation. In principle 
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the same life-time spectroscopy can be applied to heavy fermions and other exotic 
superconductors. 

The methods we proposed in this review can be initially tested by means of alter- 
native methods for determination of Cab(O), e.g. from the slope of the upper critical 
field Hc2{T) defined by the fiuctuation magnetization of the normal phase near the 
critical line, Eq. (205), being another new result derived here, or from the fluc- 
tuation conductivity of the LD model, Eq. (204). Addressing the conductivity we 
consider that the fitting to high frequency experimental data with the help of formu- 
lae (73) and (207) will give a direct method for determination of the relaxation time 
of the superconducting order parameter. A good monocrystal of layered cuprate 
or high-quality thin film with as low as possible critical temperature could ensure 
the overlap of both the suggested methods for Cooper pair life-time spectroscopy. 
At present we only know^^ that for 93 K YBa2Cu307_5 tq,* = 2to = 32 fs. We 
hope, however, that several experimental methods for determination of Ca6(0) and 
To will be mutually verified in the nearest future. Thus, the investigation of the 
Gaussian fluctuations may become a routine procedure in the materials science of 
superconductors. 

We also believe that the development of the Gaussian spectroscopy will lead 
to determination of the Ginzburg number e^j, the energy cutoff, i.e. the maximal 
kinetic energy of the Cooper pairs = cfi^ /2mab£,ab{^)- Up to now these param- 
eters of the GL theory are inaccessible. We hope that our derived self-consistent 
equation for the reduced temperature, Eq. (173), will stimulate experimentalists 
to reexamine the data for high-quality crystals in the region close to e + h o± 3% 
in order to extract e^.. Virtually all flnal results are presented by taking into ac- 
count the energy cutoff parameter c. The nonlocality corrections can be extracted 
from almost all fluctuation variables, ii e + h > 10%, but we suggest special new 
experiments to be conducted for investigation of nonlocality effects in quasi 2D 
superconductors at Tc. Analogous investigations for fluctuation diamagnetism for 
classical bulk superconductors are already classics in physics of superconductivity; 
see for example Fig. 8.5 in the well-known textbook by Tinkham.^ The universal 
scaling law for the heat capacity, Eq. (198), for the magnetization, Eq. (135), and 
conductivity, Eq. (196), versus the reduced magnetic fleld y = 2/i/c are depicted in 
Fig. 2. 

At least for conductivity the experimental confirmation for the quasi-2D super- 
conductors {h ^ r, e^,. ) can be easily achieved. How different are the animals . . . ? 
The biochemists considered that what is true for Escherichia coli holds true for the 
elephant. Analogously, we consider that B(j{Tc, B) versus B will be within 20% ac- 
curacy the same for conventional Pb layers and for strongly anisotropic underdoped 
Bi2Sr2CaCu208 in spite of the bunch of sophisticated theories of high-Tc supercon- 
ductivity. The GL theory gives the scaling law, the notions and notations, and in 
this sense the language for analysis of the fiuctuation phenomena. The precisely 
measured deviations from the GL scaling low could give the basis for further mi- 
croscopic consideration using the methods of the statistical mechanics. This is the 
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Fig. 2. Universal scaling curves ((a) semi-logarithmic plot; (b) linear plot) of a quasi-2D supercon- 
ductor for the fluctuation conductivity Ba{Tc, B) oc Ua, magnetic moment M(Tc,B) oc Um^ and 
heat capacity BC{Tc, B) oc Uq versus dimensionless magnetic field y = 2h/c. The fit of the scale 
in horizontal direction gives the GL cutoff parameter c. The scales in vertical directions are related 
correspondingly to diffusion constant of Cooper pairs / tq , effective inter-layer distance Seffi 

and 2D Ginzburg number e^., cf. Eqs. (196), (135), and (198). 



last example how the development of the Gaussian fluctuation spectroscopy could 
be of importance not only for the materials science but for the fundamental physics 
of superconductivity as well. 
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Appendix A 
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!+ Test driver program for subroutine Limes 

PROGRAM Test 



IMPLICIT NONE 

INTEGER, PARAMETER 

REAL (pr), PARAMETER 

REAL (pr) 

REAL (pr) 

REAL (pr) 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 



pr = SELECTED.REAL.KIND (30,150) 

zero = 0.0, one = 1.0 

S(0:137), C(0:137), x, xi, arg 

rLimes 

err 

N 

i 

i_Pade 
k.Pade 
is 



WRITE (* 

WRITE (* 

WRITE (* 

WRITE (* 

WRITE (* 

WRITE (* 

READ (* 

WRITE (* 

READ (* 



Test driver program for subroutine Limes 
Calculate Ln[x] 



' (15X,A) ') 
' (15X,A) ') 
' (15X,A) ') 
' (15X,A) ') 
'(A)') ' ' 

'(A)', ADVANCE='NO') ' Enter argument of Ln[x], x = ■ 

*) arg 

'(A)', ADVANCE='NO') ' Enter the number of known terms, N = ' 
*) N 



... we like this number ; -) 



IF (N > 137) N = 137 

X = arg - one 
xi = one 
is = 1 

Initialize S to store the first N+1 known partial sums 
SO, SI, 32,..., Sn-1, Sn 

Sn = X + x-2/2 - x-3/3! + x-4/4! - ... + (-D'n x-n/n! 

S(0) = zero ! *** lower bound of the subscript should start at ! *** 
DO i=l,N 

xi = xi*x 
C(i) = xi/i 

S(i) = S(i-l) + is * C(i) 



IS 

END DO 

WRITE (* 

WRITE (* 

WRITE (* 

WRITE (* 

WRITE (* 

WRITE (* 



-IS 



■(A)') ' 
'(A)') ' 



' (14X,A8,5(X,A12)) ') 'S(O) ' , 'S(l)', 

'3(4)', '3(5)' 
'(A)') ' 



'S(2)' 



'3(3)', & 



'(A)', ADVANCE='NO') 
'(F8.4,5(X,fl2.4))') 



before call: 
3(0:5) 



CALL Limes ft 

( N, ft 

3(0:N), ft 

rLimes , ft 

i.Pade , ft 

k.Pade , ft 

err ) 



call subroutine Limes to calculate Ln[x] 
in 

inout 

out 

out 

out 

out 



! Formated output of results 
] 

WRITE (*, '(A)', ADVANCE='NO') ' after call: ' 



WRITE (*, ' (F8.4,5(X,fl2.4))') 
WRITE (*, '(A)') 



3(0:5) 



WRITE (*,*) ' ' 
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WRITE (*, ' (X,A,X,F9.3,X,A,ES12.5) ') 'Ln[', arg, '] = LOG (arg) 
WRITE (*, ' (X,A,I3,A,I3,A,ES12.5) 'rLimesC, i.Pade, k.Pade, & 

'] = ' , rLimes 
WRITE (*, '(13X, lA, ES12.5)') 'err = err 
WRITE (*, '(6X, lA, ES12.5)') & 

'true error = ', ABS ( rLimes - LOG (one + x)) 



CONTAINS 



!+ Finds the limit of a series 
SUBROUTINE Limes ft 

( N, ft ! in 

S, & ! incut 

rLimes, & ! out 

i_Pade, & ! out 

k_Pade, & ! out 

err ) ! out 

Description: 

Finds the limit of a series in the case where only 
the first N+1 terms are knoun. 

Method: 

The subroutine operates by applying the epsilon-algorithm 

to the sequence of partial sums of a seris supplied on input . 

For desciption of the algorithm, please see: 

[1] T. Mishonov and E. Penev, Int. J. Mod. Phys. B 14, 3831 (2000) 
Owners: Todor Mishonov ft Evgeni Penev 



History: 
Version Date 



Comment 



1.0 



01/04/2000 Original code. T. Mishonov & E. Penev 



Code Description: 
Language : 

Software Standards: 



Declarations : 
IMPLICIT NONE 



Fortran 90. 

"European Standards for Writing and 
Documenting Exchangeable Fortran 90 Code". 



! * Subroutine arguments 
! Scalar arguments with intent (in): 
INTEGER, INTENT (IN) : : N 



width of the epsilon-table 



Array arguments with intent (inout) : 
REAL (pr) , INTENT (INQUT) :: S(0:) ! sequential row of the epsilon-table 



! Scalar arguments with intent (out) 
REAL (pr) , INTENT (OUT) : : rLimes 
INTEGER, INTENT (OUT) : : i_Pade 
INTEGER, INTENT (OUT) : : k Fade 
REAL (pr) , INTENT (OUT) : : err 

!* End of Subroutine arguments 



value of the series limes 
power of the numerator 
power of the denominator 
empirical error 



Local parameters 
REAL (pr) , PARAMETER 
REAL (pr) , PARAMETER 



! these two need no description ;-) 
zero = 0.0 
one = 1.0 



Local scalars 
REAL (pr) : : A_max 

INTEGER 
INTEGER 



maximum element of A 
index variable for columns 
index variable for rows 



! Local arrays 
REAL (pr) 



A(0:N) ! auxiliary row of the epsilon-table 



End of header 
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Parse input: the algorithm cannot employ more elements than supplied on 
input, i.e. N <= size(S) 

IF ( N > SIZE (S(:)) ) THEN 

WRITE (*, '(A)') '*** Illegal input to Limes: N > size(S)' 

STOP 1 
END IF 

Algorithm not applicable for N < 2 

IF ( N < 2 ) THEN 

WRITE (*, '(A)') '*** Illegal input to Limes: N < 2' 

STOP 2 
END IF 



I. Initialize with natural assignments 



rLimes = S(N) 

err = ABS ( S(N) - S(N-l) ) 

i_Pade = N 
k_Pade = 
A(:) = zero 
A max = zero 
k =1 



the N-th partial sum 
error -> |S(N) - S(K-l) I 
Pade appr oximant [N/ 0] 

auxiliary row initially set to zero 

max. element set to zero 

algorithm starts from the first row 



II. Main loop: fill in the epsilon table, check for convergence ... 

(provision against devision by zero employs pseudo-inverse numbers) 



DO 

IF(N-2*k+l<0) EXIT 

Update the auxiliary row A(i) of the epsilon-table 
by applying the "cross rule" . 

DO i=0, N - 2 * k + 1 

IF ( S(i+1) /= S(i) ) THEN 

A(i) = A(i+i) + one/(S(i+l) - S(i)) 
ELSE 

A(i) = A(i+1) 
END IF 
END DO 

IF(N-2*k<0) EXIT 

Update the sequential row S(i) of the epsilon-table 
by applying the "cross rule". 

DO i=0, N - 2 * k 

IF ( A(i+1) /= A(i) ) THEN 

S(i) = S(i+1) + one/(A(i+l) - A(i)) 
ELSE 

S(i) = S(i+1) 
END IF 

Check for convergence, based on A_max; see Ref . [1] 

IF ( ABS ( A(i) ) > A max ) THEN 

A_max = ABS ( A(i) ) 

rLimes = S(i) 

k_Pade = k 

i_Pade = i + k_Pade 

err = one/A max 

IF ( S(i+1) == S(i) ) RETURN 
END IF 
END DO 

k = k + 1 ! increment row index 

END DO 



END SUBROUTINE Limes 
END PROGRAM Test 
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